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TRANSPORT  PROPERTIES  OF  THE  "EXCITONIC  INSULATOR" 

# 

II.  Thermal  Conductivity 
by 

Johannes  Zittartz+ 

Department  of  Physics,  University  of  California,  San  Diego 
La  Jolla,  California 

ABSTRACT 

-  The  thermal  conductivity  of  the  excitonic  insulator'*'  is  calculated  in  the 
semimetallic  limit.  At  low  temperatures  the  main  scattering  mechanism  is  as¬ 
sumed  to  be  due  to  impurities .  Despite  a  recent  ’m  of  "supertherraal  conduc¬ 
tivity”  of  the  system  which  can  be  described  as  a  *  ondensate  of  electron-hole 
pairs  we  find  that  the  thermal  conductivity  is  ”weJl  behaved.”  It  is  shown  that 
the  thermal  conductivity  in  the  semimetallic  limit  of  the  underlying  two  band 
model  is  almost  identical  to  the  thermal  conductivity  of  a  superconductor  con¬ 
taining  magnetic  impurities  as  scatterers.  The  analogy  stems  from  the  fact 
that  the  thermodynamic  properties  in  both  cases  are  similar  as  was  discussed 
recently. 
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I.  Introduction 


In  this  paper  we  continue  the  analysis  of  transport  properties  of  the 

1  2 

"excitonic  insulator”  which  has  been  discussed  recently  in  the  literature  9  . 

In  the  preceding  paper*,  hereafter  referred  to  as  II,  we  have  calculated  the 
electrical  conductivity  of  the  excitonic  phase  in  the  semimetallic  region  where 
the  valence  bend  and  the  conduction  band  of  the  underlying  model  overlap.  The 
main  scattering  mechanism  at  low  tenjperatures  is  assumed  to  be  due  to  impurities. 
In  I  we  have  discussed  the  influence  of  the  impurities  on  the  thermodynamic 
description  of  the  excitonic  phase.  It  was  shown  that  the  situation  could  be 
described  completely  within  the  framework  of  the  Abrikosov-Gorkov  theory  origi- 
nally  developed  for  superconductors  containing  magnetic  impurities  . 

In  this  paper  we  investigate  the  thermal  conductivity  of  the  excitonic 
phase  which  can  be  described  as  a  condensate  of  bound  pairs  of  electrons  and 
holes.  The  main  reason  for  this  investigation  is -that  Kozlov  and  Maksimov  have 
recently  claimed  "superthermal  conductivity"  for  such  a  system^.  Such  a  novel 
situation  certainly  would  be  of  major  experimental  interest.  However,  we  shall 
show  in  the  course  of  this  work  that  the  thermal  conductivity  of  the  excitonic 
phase  is  "well  behaved,"  in  fact  it  is  similar  to  the  thermal  conductivity  of 
a  superconductor  containing  paramagnetic  impurities  as  scatterers^.  The  contra¬ 
dictory  claim  of  reference  5  is  based  on  the  Landau  criterion  for  superfluidity; 
the  derivation  of  the  criterion  for  the  excitonic  phase  contains  a  serious 

7 

mathematical  error'. 

To  start  with  we  shall  first  discuss  Kubo's  formulas  for  the  transport 

coefficients.  Measuring  energies  relative  to  the  chemical  potential  the  elec- 

•  8 

trie  current  j  and  the  heat  current  U  are  defined  as  follows 


j  -  -f-(;r-£vT)  -£vt 

s  '-¥<Zr--$r)~ 7?t 


where is  the  chemical  potential,  T  the  temperature  and  e  the  electronic 

i 

charge;  <S  is  the  d.c.  conductivity.  The  thermal  conductivity  8 fc  is  obtained 
from  (l)  by 

U,  -  -“it  VT  ;  j  cO  (2) 

which  leads  to 


The  d.c.  conductivity  S'  has  been  calculated  in  the  preceding  paper  whereas  the 

o 

transport  coefficients  and  L,  are  given  by  (we  anticipate  spherical  syn- 
metry) 

oo  p 

L  *  J  U*  <  U(o)  U(t+.A)> 

,  ~  0  <*> 

L,  sxJcfiJolA<^(o)^;A)>  . 

&  0  • ««*> 

Here  U  and  n  are  the  Heisenberg  operators  for  the  energy  and  electric  currents, 
respectively.  Hie  brackets  indicate  a  thermal  as  well  as  an  average  over  im¬ 
purities  as  discussed  In  I  and  II. 

The  two  band  model  we  shall  use  has  been  discussed  in  detail  in  I  and 
the  preceding  paper  (II).  Hie  single  particle  energies  in  the  two  bands  are 
assumed  to  be  given  near  the  Fermi  surface  by  (see  eq.  (II. 2) 

Vp>--Vp)»-U  *  (5) 
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vber*  for  mathematical  convenience  ve  assume  equal  band  masses.  We  shall  see 


later  on  that  In  the  semimetallic  limit  (large  pQ)  the  thermoelectric  coef¬ 
ficient  Lf  is  smaller  by  a  factor  WjL  as  compared  to  L  .  As  we 

/  Zm 

shall  systematically  neglect  such  terms,  we  confine  our  attention  to  the  cal¬ 
culation  of  the  coefficient  L  . 

Another  useful  representation  of  L  can  be  derived  from  the  Kubo  form¬ 
ula  (4)^.  Consider  the  causal  correlation  function 

* -jr  <T  ui  v  u«o>  <«> 

where  ^  are  imaginary  times  with  0  i  T,  fi  •  Fourier  transforming  ac¬ 
cording  tc 

-  f  Z  WO  e'*™  f  v  siniTp"1  c, 

we  continue  the  discrete  imaginary  frequencies  IV  to  the  full  complex  energy 
plane.  Ihe  coefficient  is  then  given  by 

L  »  £  «»»  i  [T(«V,f)-? (»-;S)l  *  tiw,  Ovn  1W«  J)  (®) 

W^>0  t0-%0  w 

where  4  is  a  positive  infinitesimal. 

However,  we  have  to  be  careful  in  using  this  formula  for  answering  the 
question  of  "superthermal  conductivity.”  This  is  because  the  mixing  of  "retarded 
and  "advanced"  response  in  formula  (8)  neglects  the  possibility  of  an  anomalous 
behavior  of  the  real  part  of  in  the  limit  •  Similar 

to  the  problem  of  infinite  conductivity  in  superconductors,  we  expect  that  the 
possibility  of  an  undamped  heat  current  which  could  be  accelerated  by  the  tem¬ 
perature  gradient  would  show  up  in  the  real  part  of  the  retarded  correlation 
function  T(u>-ti£l  (thus  leading  to:  ~  U,  *■  COhst,  ).  Instead  of 

'  'XJ  *“T 


formula  (8)  we  shall  consider  the  following  expression 


L 


?(<*+;$)  -  c 

\  UJ 


(9) 


where  the  constant  C  will  be  determined  later  on  by  comparison  with  the  cor¬ 
responding  expression  for  the  normal  state  behavior.  Formula  (9)  bears  resem¬ 
blance  to  the  familiar  expression  for  the  electrical  conductivity  (see  preced¬ 
ing  paper).  ?(<*♦! would  correspond  to  the  paramagnetic  contribution  to  the 

p 

electric  current  (denoted  by  K  in  II)  whereas  the  constant  C  would 

be  equivalent  to  the  diamagnetic  contribution  (denoted  by  in  II)10. 

In  section  II  we  derive  a  tractable  form  for  not  dwelling  on 

additive  constants  which  are  assumed  to  be  absorbed  in  the  constant  C.  In 
section  III  we  make  an  ansatz  for  the  constant  C  based  on  the  required 
form  in  the  normal  state.  We  also  shall  be  guided  by  the  corresponding  expres¬ 
sion  for  the  diamagnetic  contribution  in  the  problem  of  electrical  conductivity. 
The  resulting  expression  for  the  thermal  conductivity  will  turn  out  to  be  "well 
behaved"  as  U) 0  .  Thus  there  is  no  "superthermal  conductivity."  In  section 
IV  we  calculate  the  final  expression  in  analogy  to  calculations  in  the  preceding 
paper.  Section  V  contains  explicit  evaluations  in  several  limits. 
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II.  Derivation  of  the  heat  current  correlation  function 


The  model  has  been  discussed  in  detail  in  papers  I  and  II  (see  section 
II  of  preceding  paper).  Denoting  creation  and  annihilation  operators  for  a- 
and  b-  electrons  by  U+  and  respectively  we  have  introduced  a 


Itaribu  notation 


▼S?  -( -)  /  w  -(!>;»;)  . 

\  ftp/  MM 

The  electric  current  operator  has  been  written  as  (eq.  11.19) 

j(t)  *  el  y(pz)  *<p)  TKpt) 

A  '  * 

MM,  r  HM 


where  the  velocity  matrix  'tf(p)  is  explicitly 

MM 

....  ( y?  M  .  - . 


*(?)  * 


0  V 


)  =  if  V£  . 


Similarly  the  spatially  uniform  heat  current  operator  can  be  defined  by 


U(t)  (13) 

•—  £  "*■  "*  - 

where  are  imaginary  times  and  £  a  positive  infinitesimal  indicating 

the  ordering  of  operators.  Expression  (13)  is  the  appropriate  form  when  deal¬ 
ing  with  imaginary  times  in  which  case  the  Heisenberg  field  operators  are  given 
by  equation  (II. 8)  of  the  preceding  paper.  Introducing  (13)  into  the  correla¬ 
tion  function  (6)  we  obtain 
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(±.. 

-M 

/_L 

( 

*  v(p<)  z(p>  v?y  yW  z(i!  yw>  I 

I  ** 


a*) 


\  =  Vy 


Actually  the  time  derivatives  in  this  expression  should  have  been  included  in 
the  brackets,  as  they  do  not  commute  with  the  time  ordering  operator  T.  By 
taking  these  derivatives  outside  we  have  neglected  a  term  proportioned  to  the 
S  -  function  .  When  Fourier  transformed  this  correction  term 

gives  a  constant  which  we  assume  to  have  been  included  in  the  constant  C  in 
equation  (9)» 

The  further  calculation  proceeds  in  a  way  similar  to  the  calculation 
of  the  electrical  current  correlation  function  in  the  preceding  paper.  In¬ 
troducing  Gresns  functions  (see  I  and  II)  and  neglecting  vertex  corrections 
for  a  moment  we  get 


P  ( v\) 


_  3L  )/jL  .  IL.  \ 

•a.  I  c>T dzj  ’ 

Z  Trace  v(p)  fy)  I 

p  *+  «M.  »  I  • 

-  'x.-va 
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Fourier  transforming  we  obtain 


?(■'?)  m~'h  IT"!"*  2LTmce  'tty  G(MatT**)  &  ^  ^l6) 

^  d“-o 

where  the  time  derivatives  have  been  written  as  derivatives  with  respect  to  j  , 
Coopering  this  expression  with  the  electrical  current  correlation  function 
K?(iv)  (11,21)  we  see  that  both  are  identical  except  for  the  prefactor, 
the  J  -exponential  and  g  -derivatives  and  except  that  we  have  left  out  the 
vertex  correction  for  the  velocity  matrix  in  (16).  The  vertex  corrected  ex¬ 
pression  T/Jj^may  be  written  down  immediately  from  this  comparison  with  -(11,21) 


where  W  has  been  calculated  in  section  III  of  the  preceding  paper.  We  remark 
that  as  usual  the  U>  -summation  has  to  be  performed  first  as  the  double  sum¬ 
mation  converges  only  conditionally.  Expression  (IT)  is  our  final  form  of  the 
heat  current  correlation  function.  The  constant  C  which  according  to  equa¬ 
tion  (9)  has  to  be  subtracted  from  (17)  will  be  considered  in  the  following 
section. 


III.  Determination  of  the  subtraction  constant 

We  shall  determine  the  subtraction  constant  C  in  formula  (9)  by  anal¬ 
ogy  from  the  corresponding  constant  in  the  electrical  conductivity  problem, 
i.e.  the  diamagnetic  contribution.  The  latter  is  given  from  equation  (II.18) 
of  the  preceding  paper 


IjVpO  vr> +  (Y,2vp>)  \<r>  ] 


where  the 


t  ^are  the  general  band  energies,  Hftfc(p)  the  occupation  numbers  in 

'  i  ^ 


-8- 


momentum  space.  We  rewrite  the  occupation  numbers  in  terms  of  the  spectral 
functions*2 


(19) 


where  is  the  Fermi  distribution  function^ and  the  spectral  functions  are 
expressed  through  the  one-particle  Greens  functions 


Sft/p.w)  =  '  ,»*•$)  -Q/pl'*>-''3)]  (2°) 

Thus  we  can  write  for  (18) 

-K*-ir f e>)54 </>.“)  +(‘”")  .  <a> 

The  product  'f(to)S(p|w)  can  be  Interpreted  as  the  density  of  particles  with  no- 

M* 

mentum  p  and  energy  <4  . 

It  is  intuitively  clear  that  in  the  thermal  conductivity  problem  we 
should  replace  the  charge  t  in  expression  (21)  by  the  energy  co  .  Therefore 
we  make  the  following  ansatz  for  the  constant  C: 

c“  £ -a-.*) 

»-  (22) 

=  V ?  IF  1  I  +  (b-> 

o  r  -  e  +1  -  ^  U 

where  in  the  second  line  we  have  explicitly  introduced  the  functions  and  5  . 
The  (o  -factor  is  generated  by  the  ^-derivatives,  and  ^  is  a  small  positive 
number  which  has  to  be  set  equal  to  zero  finally  as  indicated  in  (22). 

The  -integral  can  be  transformed  into  a  simulation  over  imaginary  fre¬ 
quencies  1(4^  .  This  is  done  by  closing  the  contour  of  the  integral  for  the 

first  integrand  in  (22)  by  a  large  semicircle  in  the  upper  half  plane  and  for 
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the  second  integrand  in  the  lower  half  plane.  Siamring  over  the  residues  at  the 
poles  of  the  Fermi  function  at  we  get  at  once: 

C  *<>~)  <*> 
£  m  **  ^  J 

where  the  <jo -summation  has  to  be  performed  first  as  usually.  In  order  to  make 
contact  with  the  correlation  function  expression  (IT),  the  right  hand  side  in 

(23)  is  further  transformed  by  a  partial  p  -  integration  over  the  Brillouin 

m, 

zone.  In  analogy  to  the  calculation  in  the  preceding  paper  (eq.  11,22)  we  get 


••  1+  ** 


ar=o 


(24) 


where  are  the  group  velocities  ~  (12)*  Using  the 

Greens  functions  of  the  preceding  paper  (11,9),  we  have  after  taking  the  deri¬ 
vatives  and  combining  al^  terms  in  analogy  with  expression  (II. 23): 


For  better  comparison  with  ?(.V)  (IT)  we  also  have  substituted  -+  iff 
Comparing  this  expression  with  expression  (IT)  for  the  correlation 
function  ?(oj)  we  see  that 


C  =  ?(i»*o) 


(26) 


because  the  vertex  corrected  velocity  matrix  W  in  (IT)  reduces  to  the  free 

A  ** 

vertex  function  V  at  VrO  (see  equations  (II. 25, 26,30)  of  the  preceding 

*• 

paper).  As  the  statement  (26)  is  independent  of  the  fact  whether  the  system 
is  in  the  normal  state  or  the  excitonic  state,  we  can  conclude  quite  generally 
that  the  thermal  conductivity  is  "well  behaved."  This  is  because  (26)  guarantees 
that  the  real  part  of  the  numerator  in  (9)  vanishes  in  the  limit  U)-*0  .  This 
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conclusion  will  be  confirmed  by  the  calculations  in  the  next  two  sections. 


IV.  Evaluation  of  the  thermal  conductivity 

Combining  the  two  expressions  (17)  and  (25)  we  can  follow  very  closely 
the  calculations  of  the  corresponding  expressions  ( II, 21,23)  in  the  preceding 
paper.  We  first  perform  the  momentum  integral  evaluating  the  toted,  expression 
at  the  Fermi  surface.  In  analogy  to  the  result  (II, 33)  of  the  preceding  paper 
we  get 


p(o>)  *?M-C  = 


L  It  tr  •  •  a 

""v’) 


(27) 


where  jp  is  the  density  of  electrons  or  holes,  respectively,  and  the  function 
F  has  been  defined  in  equation  (11,34).  The  further  calculation  differs  slightly 
from  the  corresponding  calculation  in  section  IV  of  the  preceding  paper  because 
of  the  presence  of  the  ft -exponential  in  (27).  We  write  the  frequency  summation 
as  a  contour  integral  in  the  complex  plane: 


>)jf 


Is  -  -jfe  J**  V-  F(z, *-'»>)  f 


(28) 


the  path  Cq  being  defined  in  equation  (11,35)  of  the  preceding  paper.  Because 
of  the  y  -exponential  we  have  to  choose  the  Fermi  function  in  (28)  instead  of 
the  tanh- function  in  (11,35)  in  order  to  insure  the  vanisning  of  the  integrand 
for  large  z.  Ihe  integral  (28)  is  further  transformed  into  an  integral  over  the 
discontinuities  of  the  function  Ffe,  2-iy)  across  its  cuts.  Ihese  cuts  are 
defined  in  equation  (11,37).  The  contribution  f  rom  the  first  cut  is  given 
by 


2 


u>.  6+1  J 


(29) 
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where  wi  is  the  gap  in  the  excitation  spectrum  of  the  system  (II.A10).  Simi- 

o 


larly  the  contribution  from  the  cut  gives 
,  'v-“e 

*4  2L  J  0,1  gjj£  )~F(l,Z-i9-\S)l 

iy-*o  er  +  I  -* 


er  ♦  l 

-(Zx-iy)^ 


2  ->  iy-x 


»„  ft  + 1 

where  we  have  used  (II.  38).  Ihe  contributions  from  the  two  ^uts  Cg  and 
can  he  transformed  to  integrals  over  the  cut  in  the  same  way.  Combining 

all  4  terms  and  rearranging  we  obtain  for  the  right  hand  side  of  (27) 

gL  00 

p(,v)  * '  w  v*/-x{CRW **)tj 

I’lf’fx+I  <5J  **  1 »)  “ X + f v)]£ta H  h  C0Sh/2*ii)>)fr  I 

According  to  formula  (9)  we  have  to  take  the  analytical  continuation  of 
this  expression  from  *»Jy  to  2*u)+id  and  then  pass  to  the  limit  U>-*0.  IMs 
process  leads  to  the  following  expression  for  the  transport  coefficient  1 

00 

L  *  zfa  Jc**^F(xti3/)i-u>-i5)^anVi'fc  coih(2x-«3)j  -s<yih(ZK-ui)fiJ 

0  (a)^ 


w-/ 

■Rx-i?;xnoti^tanlr^  005(1(2^0))^  -si*J)(2x+u>)yj  J'  j 


where  we  have  left  out  the  contributions  from  F(x-iffx-w-i  f)  and 
which  according  to  (11.42)  vanish  in  the  limit  U>~*0  .  Furthermore  we  can 

transform  in  the  first  part  of  the  integral  and  split  the  integration 

into  two  terms.  The  first  of  these  terms  is  the  integral  from  (0-10  to  U)^ 
and  gives  for  small  (a) 
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«0 

co”  J  R^iMi  jjtofnh  cosh(ixtv))^  -  sinh(2xn»)^ 
'V  FfWet'i^  t*0-i CO$)](2uft^)  -  sitih(Z<Oojf)J  —  0 


as  discussed  in  equations  (11.44,45)  of  the  preceding  paper.  The  second  inte¬ 
gral  (<0e£x<«o)  can  be  combined  with  the  second  part  in  (32)  leading  to 

a*-  ©o  . 

L  *  £  >3  s  *>  |  “*<*5)  /  »> 

where  we  have  replaced  . 

Comparing  this  expression  with  equations  (11.46,47)  of  the  preceding 
paper  we  see  that  the  d.c.  conductivity  apart  from  the  factor  &  may  be  ob¬ 
tained  from  (34)  by  leaving  out  the  y -derivatives.  This  is  explained  by  the 
fact  that  the  heat  current  (13)  differs  from  the  electric  current  (11)  by  the 
time  derivatives  which  eventually  are  ce presented  by  the  y  -derivatives  in  the 
correlation  function  expression  (IT) •  In  the  same  way  it  is  clear  that  the 
thermoelectric  coefficient  Lt  (  4)  is  given  by  (34)  with  one  ^-derivative  only, 
as  there  is  only  one  heat  current  involved  in  the  correlation  function  (4)  for 

L  13 

Therefore  obviously  vanishes.  The  reason  for  this  is  that  we  have 
evaluated  the  momentum  integral  at  the  Fermi  surface  in  passing  from  expression 


(17)  to  expression  (27).  Correction  terms  would  be  smaller  by  a  factor 


K.TM 


which  we  have  neglected  systematically.  Thus  the  thermoelectric  coefficient 
L,  is  zero  to  a  sufficient  accuracy. 

The  thermal  conductivity  at  is  therefore  give;i  from  (3)  and  (34) 
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at  Jdx  x1**^ 


hfx) 


-rj  +2p'o-h) +rtr  (35) 

where  we  have  Introduced  the. function  ^  F(HtJJjK-if)  from  equation  (llA8) 
of  the  preceding  paper;  the  different  terms  and  symbols  are  also  explained  in 
section  V  and  the  Appendix  of  the  preceding  paper.  Explicit  calculations  in 
several  limits  are  considered  in  the  final  section. 

V.  Discussion  and  explicit  calculations 

We  may  compare  our  final  expression  (35)  with  the  corresponding  formula 
for  the  thermal  conductivity  of  a  superconductor  containing  magnetic  impurities 
as  derived  in  reference  6  (equation  (3.1)^  The  two  expressions  are  identical  if 
we  replace  the  total  density  of  particles  j)  in  the  formula  of  reference  6  by 

in  our  case  which  is  due  to  the  two  types  of  carriers.  The  similarity 
is  not  surprising  as  we  have  seen  in  reference  2  (paper  I)  that  the  excitonic 
phase  in  the  presence  of  impurities  has  the  same  thermodynamic  properties  as  a 
superconductor  with  magnetic  impurities.  Therefore  one  expects  that  the  reac¬ 
tion  to  a  temperature  gradient  is  similar  in  both  systems  because  the  different 
charges  of  the  bound  pairs  are  not  significant  in  this  case. 

The  fact  that  the  expression  (35)  is  well  behaved  over  the  whole  tem¬ 
perature  range  clearly  disproves  the  claim  of  reference  5  that  the  excitonic 
phase  should  have  "superthermal"  properties.  It  remains  to  discuss  the  behavior 
of  the  thermal  conductivity  in  several  limits  which  are  accessible  to  explicit 
calculations . 

In  the  normal  state  (A=0)  we  obtain  from  (35)  using  (II.51)  of  the 


preceding  paper 


where  =.  P 

ir.  ir. 

parison  with  the  "normal" 


is  the  transport  collision  time  (II. 50).  Upon  corn- 
conductivity  (11.52)  we  have  the  Wiedemann-Franz  law 


(37) 


Using  (36)  we  write 
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3  £j\r. 

2.1T2- 


J  dx  xlsethl& 
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_ hU) _ 

z[^iZ?-r]+zr'(i-h)4rir 


(38) 


In  order  to  evaluate  this  expression  asymptotically  near  T  =  0  we  have  to  dis¬ 
tinguish  between  the  two  cases  whether  the  excitation  spectrum  of  the  system 
has  a  gap  ( o)  or  not.  In  the  first  case,  i.e.  for  low  inpurity  con¬ 
centration  (  o(  =  ^  <  |  ,  see  the  preceding  paper),  we  obtain,  using  the 

A 

expressions  (11.54)  and  (11,55)  of  the  preceding  paper,  respectively 


T  >0 


(39) 


Uius  the  thermal  conductivity  goes  to  zero  exponentially  as  long  as  there  is  a 
gap  in  the  excitation  spectrum. 

In  the  gapless  region  (*>0  we  get  at  T  =  0  using  (11,58)  of  the  pre¬ 
ceding  paper  and  starting  the  integration  in  (38)  at  U)0«0 


-x 


-2- 


,  *r' 

,+7T« 

ir. 


(4o) 


Comparing  with  (II. 59)  we  see  that  in  the  gapless  region  the  ratio  ^7-*  is 

identical  to  the  ratio  for  the  d.c.  conductivity,  therefore  the  Wiede- 

►t 

mann- Franz  law  holds,  too. 
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Finally,  we  consider  the  transition  temperature  region  where  the  system 
Is  always  gapless  /  <U  A o)  •  Expanding  (38)  In  powers  of  A. 

we  conclude  In  analogy  to  equations  (II. 62,65)  of  the  preceding  paper  that 

7  -I  -  const  t=j<  j  W 

/  «c 

Huserical  solutions  for  the  superconductor  case  In  the  whole  temperature  range 
are  reported  In  reference  6. 
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